
AMENDED SPECIFICATION (marked-up) 

METHOD AND ELECTRONIC DEVICE USED TO SYNTHESISE THE SOUND OF 
CHURCH ORGAN FLUE PIPES. BY TAKING ADVANTAGE OF THE PHYSICAL 
'DELING TECHNIQUE OF ACOUSTIC INSTRUMENTS. 



applicable. 



Not applicable. 



Not applicable. 



Not applicable. 



CROSS-REFERENCE TO RELATED APPLICATIONS 



STATEMENT REGARDING FEDERALLY SPONSERED 
RESEARCH OR DEVELOPMENT 

THE NAMES OF THE PARTIES 
TO A JOINT RESEARCH AGREEMENT 



INCORPORATION-BY-REFERENCE OF MATERIAL 
SUBMITTED ON A COMPACT DISC 



BACKGROUND OF THE INVENTION 

( 1 ) Field of the Invention 

The present patent application refers to a method and electronic device used to synthesise the sound of 
church organ flue pipes, by taking advantage of the physical modeling technique of acoustic 
instruments. 

Numerous numerical algorithms of physical-mathematical models have been developed based on the 
examination of the physical behaviour of organ flue pipes and the sound they produce, in order to 
synthesise the sound emission of aerophone instruments in real time. These models are based on the 
mutual symbiotic interaction between a non-linear active section, generally defined as "excitation", 
and a linear passive section, generally defined as "resonator". The relative numerical algorithm 
extemporarily produces a sequence that represents the sound of the instrument analysed and translated 
into a physical model. The sound is characterised by an initial time interval, defined as "attack 
transient", during which intensity increases up to a certain value. The intensity value is indefinitely 
maintained over time during the second phase, defined as "sustain phase", during which the waveform 
is approximately periodic. The analytical characteristics of this waveform, of which the most 
important is fundamental frequency, depend on each of the parameters that regulate the operation of 
the numerical simulation. 

Being the simulation performed in the time domain instead of the frequency domain because of the 
presence of numerous non-linear functional blocks, the relation between the set of parameters and 
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each spectral characteristic of the generated sequence is extremely difficult to establish a priori. 
The characteristics can be altered by changing the set of parameters, often empirically, and then 
evaluating the effect of such a change a posteriori. In particular, the fundamental frequency also 
depends on the quantitative characteristics of excitation, and not only on the frequency response of the 
resonator; being the evolution of the sequence extremely chaotic during the attack transient phase, the 
phase of the fundamental frequency cannot be pre-determined once the sustain phase has been 
reached. These two peculiarities are unacceptable in high-polyphony electronic musical instruments, 
such as church organs. 

BRIEF SUMMARY OF THE INVENTION 
The present invention consists in an audio-digital synthesis system based on digital signal processors, 
which contains a programme of physical simulation of the sound generation of organ flue pipes. The 
programme is divided into three fundamental, conceptually independent sections: the first section 
generates the harmonic part of the sound; the second section generates the aleatory part of the sound; 
the third section processes these components by means of a transfer function with two inputs and one 
output, thus obtaining the sequence that represents the sound of the organ pipe. Because of the 
independence of the section that generates the harmonic part of the sound, the fundamental frequency 
and the phase of the whole waveform generated by the programme can be determined a priori. 
The numerical parameters of the simulation programme are partially contained in a static memory and 
partially obtained by processing information from an electronic musical keyboard and from a set of 
user controls in real time. They determine the fundamental characteristics of the generated sound, 
among which the main characteristics are pitch, intensity, time envelope, harmonic composition and 
aleatory component. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS 
OF THE DRAWINGS 

For major clarity the description of the method and device according to the present invention continues 
with reference to the enclosed drawings, which are intended for purposes of illustration only and not in 
a limiting sense, whereby: 

- Figure 1 shows a realisation of a digital electronic musical instrument used to synthesise sounds of 
musical instruments by taking advantage of the physical modeling technique of the invention. 
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- Figure 2 shows the three fundamental functional blocks and relative interconnections of an audio 
digital synthesis programme of the sounds of church organ flue pipes according to the invention. 

- Figure 3 shows a flow chart that explains one of the three blocks of Fig. 2, according to which a 
sequence that represents the harmonic part of the sounds of church organ flue pipes according to the 
invention is generated. 

- Figure 4 shows a stable realisation of a digital harmonic oscillator with two status variables 
according to the invention. 

- Figure 5 shows a procedure used to generate the time variation of the operational frequency of the 
harmonic oscillator shown in Fig. 4 according to the invention. 

- Figure 6 shows a flow chart used to generate the aleatory component of the time progression of the 
operational frequency of the harmonic oscillator shown in Fig. 4 according to the invention. 

- Figure 7 shows an example of time envelope used in the generation of the sequence that represents 
the harmonic part of the sounds of flue pipes according to the invention. 

- Figure 8 shows a flow chart of a low frequency oscillator used in the generation of the sequence that 
represents the harmonic part of the sounds of flue pipes according to the invention. 

- Figure 9 shows a time progression composed of non-rectilinear sections, according to which the 
frequency of an oscillator can be changed without perceiving an alteration of timbre pitch according to 
the invention. 

- Figure 10 shows an algorithm for the generation of a pseudoimpulsive periodic sequence according 
to the invention. 

- Figure 1 1 shows a set of interconnected functional blocks that explains one of the three blocks of Fig. 
2, according to which a sequence that represents the aleatory part of the sounds of church organ flue 
pipes according to the invention is generated. 

- Figure 12 shows a status device used to limit the difference between two consecutive samples of a 
sequence according to the invention. 

- Figures 13 and 16 show an example of wave envelope used during the attack transient phase of the 
generation of sounds of flue pipes according to the invention. 

- Figure 14 shows a wave envelope used to generate the aleatory component of the sounds of flue pipes 
according to the invention. 
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- Figure 15 shows an architecture that explains one of the three blocks of Fig. 2, representing a 
mathematical model of the resonator of the church organ flue pipes according to the invention. 

- Figure 1 7 shows the mutual interaction between two functional blocks necessary for the realisation of 
a generic harmonic oscillator according to the invention. 

- Figure 18 shows an example of a pseudoimpulsive periodic waveform generated by the algorithm of 
Fig. 1 0, used to generate the aleatory component of the sounds of flue pipes according to the invention. 

- Figure 19 explains the operation of the status machine of Fig. 12 according to the invention. 

DETAILED DESCRIPTION OF THE INVENTION 
With reference to the aforementioned figures, the electronic musical instrument of the invention is 
physically composed of a set of components, whose type, arrangement and interconnection are shown 
in Fig. 1. 

The embodiment is shown for mere illustrative purposes, since it neither represents the central 
innovation element of the present patent nor the only and necessary realisation of an electronic musical 
instrument used to synthesise the sound the organ pipes by means of algorithms of 
physical-mathematical simulation. With reference to Fig. 1, the information from a musical keyboard 
(1) and a set of user controls (2) is processed by a control unit (3), which regulates the operation of a 
DSP (6) by means of a plurality of numerical parameters contained in a ROM (4). The DSP (6) 
executes the synthesis programme of the sound of the organ pipe in real time, upon management from 
the control unit (3), using a RAM (5) to write and read temporary data. The product of the synthesis 
programme is a numerical sequence that is suitably converted by a DAC (7) into the analogue signal 
representing the sound of the organ pipe, which can be reproduced with an amplification system and a 
loudspeaker (8). The synthesis programme, which is the central innovation element of the present 
patent, includes three sections. Each section has a fundamental function in the numerical simulation of 
the sound emission of the organ pipe, as shown in Fig. 2. 

The block (9) generates a main harmonic sequence (10) composed of a set of harmonic lines, whose 
amplitude and frequency conveniently change over time. By using this sequence and taking advantage 
of part of its composition, the block (11) generates a pseudoaleatory signal that represents the chaotic 
component of the sound. The aforementioned sequences are the two input signals of the linear 
resonator (12) that models the frequency response of the resonant part of the multiple qualities of 



4 



Serial No. 10/549,388 
Docket No. 024783.003 

organ flue pipes, and whose output (13) is the sequence that represents the sound of the organ pipes. 
The block diagram of Fig. 3 is a detailed view of the functional blocks of the harmonic component 
generator (9). The oscillator (14) generates an approximately sinusoidal waveform (16). The 
fundamental frequency of the waveform changes over time within a range of values comprising the 
fundamental frequency of the generated musical note. The details of the embodiment of the oscillator 
and the criterion used to change frequency over time are illustrated below. 

The waveform (17) is obtained from the sequence (16) through the non-linear block (15): if the 
sequence (16) were exactly a sinusoidal sequence 

x[n] = sin[cQ 0 n], 

the sequence (17) would be 

y[n] = 2 (sin[co 0 n]) 2 - 1 = -cos[2co 0 n] = sin[2co 0 n - n/2] 9 
that is to say a sinusoid with double frequency than the sequence (16). 

Each of the two sequences (16) and (17) is amplified by the relevant multipliers (18a) and (18b), and 
limited by the functional blocks (19a) and (19b) to values within the ±CLIP1 and ±CLIP2 intervals. 
The outputs of the blocks (19a) and (19b) are multiplied by two sequences produced by the envelope 
generators (20a) and (20b), respectively, as illustrated below, and the resulting products are summed 
to the node (21). The sum is a sequence produced by a series of linear and non-linear instantaneous 
operations performed on the waveform (16). If the waveform were exactly the sequence x[n], a 
sequence would be obtained in the node (21) whose spectrum would be formed by harmonic 
components multiple of © 0 (including co 0 ). 

As illustrated below, the sequence (23) is a low frequency waveform, whose purpose is the amplitude 
modulation of the harmonic sequence through the product (22). 

The element (24) is a delay line whose impulse response is the sequence 5_i[n - N], Together 

with the products and the sum of the block (25), this element forms a linear filter whose impulse 
response is 

CBYP + CDELz" N . 

The block (26) is a non-linear instantaneous function described by the following formula: 

Xx) = (x + xo)-(x + xo) 3 + y 0 
where xq and y 0 are independent parameters. The purpose of the block is to modify the mutual 
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proportion between the amplitudes of the harmonic components of the sequence processed by the 
block. 

The block (27) is a band-pass filter, whose peak frequency corresponds to the fundamental frequency 
of the input sequence. The parameter Q of the filter is tuned up to obtain the only fundamental 
frequency of the input harmonic sequence with excellent approximation. Moreover, being the phase 
response of the filter null in correspondence of the peak frequency, the phases of the fundamental 
frequency of the input and output signals of the filter are equal. This characteristic enables to sum the 
input and output sequences of the filter, with no elision effect in the fundamental frequency: the block 
(28) sums the sequences (weighing them with the parameters GAIND and GAINF), in order to alter 
the proportion in amplitude between the fundamental harmonic component and the group of all other 
harmonic components. The output of the block (28) is the main harmonic sequence (10). 
The sinusoidal oscillator (14) consists in a special embodiment of the ordinary harmonic oscillator 
with two status variables, with necessary measures to improve the robustness to the variation of the 
operational frequency in real time. 

Fig. 4 shows the cycle of operations performed at each sampling interval on the two conveniently 
initialised status variables VAR1 and VAR2. The parameter F determines the frequency of the 
sinusoid produced by the status variables oscillator that is composed of the steps (29) and (31) in the 
ordinary configuration. The disadvantage of the ordinary configuration is that it cannot suffer 
variations of the parameter F in real time without altering the amplitude of the sinusoids described by 
the same variables, in function of the current value of the status variables. Moreover, depending on the 
numerical precision of the oscillator's status variables, reductions of the oscillation amplitude can 
occur even in stationary conditions. It is sufficient to amplify the variable VAR2 by a factor 1+e (with 
e positive, but close to zero) by means of the step (30) and limit the width of the variable VAR1 by 
means of the step (32) to values within the interval ±1. Using these measures, the variable VAR1 
describes a unitary amplitude sinusoid with excellent approximation. This variable is the output (16) 
of the block (14) of Fig. 3. The parameter F depends on the frequency f according to the relation 

F(f) = 2sin(7if7f sr ) 

where f sr is the sampling frequency. The frequency f can vary in real time within an interval [f 0 - Af, f 0 
+ Af] sufficient to have the frequency changes perceived, without a collateral amplitude alteration. 
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Having defined the deviation from the central frequency fo as 8f, this parameter changes in real time 
according to the scheme of Fig. 5. Likewise the signal (23), the signal (33) is a low frequency 
waveform whose purpose is the frequency modulation of the generated sinusoid; with the support of 
the variable VAR1, the block (34) generates an aleatory waveform of "sample and hold" type, 
according to the scheme of Fig. 6. Ultimately, 6f varies according to a constant PITCH parameter 
(which, assuming a value in an arbitrary interval [1-5, 1+8], determines the fine tuning of the sinusoid) 
of an oscillating sequence (33) and the aleatory sequence (34). The block (34) is described in Fig. 6: 
every time the variable VAR1 passes from a negative value to a positive value, the variable 
RNDPTCH is updated to a new value NEWRND, which is an aleatory variable with a probability 
density function uniformly distributed in the interval [l-SRNDP, 1+8RNDP], being 8RNDP an 
independent parameter. 

The two generators (20A), (20B) produce two 5-segment envelope signals, whose progression is 
generically illustrated in Fig. 7. Tl . . .T4 are the time intervals in which the signal passes from level L0 
to LI, from LI to L2, from L2 to L3 and from L3 to zero, respectively. The generators start 
producing the respective envelope signals upon a "note on" event. Level L2 is maintained over time 
for an indefinite interval SUSTAIN, whose end coincides with the corresponding "note off' event. 
Each of the two generators uses its own set of these 8 parameters. 

The signals (23) and (33) are produced by a "Low Frequency Oscillator" shown in Fig. 8. The 
generation method of the triangular waveform with unitary amplitude and frequency TRFREQ 
illustrated in the block (35) is implicit. The parameters TRFREQ, TRAMPL, TROFFSET, 
TRCOEFF1 and TRCOEFF2 determine the conformation of the two signals (23) and (33), whose 
common fundamental frequency is TRFREQ. In particular, the signal (32) is a triangular wave of 
average value TROFFSET and semi-amplitude TRAMPL, while the signal (33) is formed by sections 
of parabolas, as shown in Fig. 9. The relation between the values TRCOEFF1, TRCOEFF2 and the 
independent parameter K is biunique. The special progression of the signal (33) is necessary to obtain 
a triangular frequency modulation as exactly as possible (ref. Fig. 5) around the nominal frequency fo, 
with equal progressions of the positive and negative semi-periods, if they are expressed in semitone 
cents. 

The architecture of the generator (1 1) of Fig. 2 is illustrated with details in Figures 10, 1 1 and 12. With 
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reference to Figs. 3 and 10, the signal (16) produced by the sinusoidal oscillator (14) is amplified by a 
factor RTINGAIN, limited in amplitude by the block (36) to values within the interval ±1, and then 
processed by the high-pass filter (37). Finally the non-linear block (38) cuts the signal's negative 
values. At the output of the block (38) the signal (illustrated in Fig. 13) produced by the envelope 
generator (39) is summed and the result is multiplied by the parameter RTGAIN. The result RATE is a 
sequence of values used in the non-linear block (42) defined as "RATE LIMITER", which is part of 
the structure described in Fig. 11. With reference to Fig. 11, the functional block (40) generates a 
white aleatory sequence, with a uniformly distributed probability density function processed by the 
low-pass filter (41). The obtained sequence is the input signal of the structure formed by the delay 
lines NBDL1, NBDL2, NBDL3, NBDL4, the sums NBS1, NBS2, NBS3, the multipliers NCGAIN, 
NBFBK and the non-linear block (42). The set formed by these elements, including the topology of 
interconnections, is defined as "NOISE BOX". The signal generated by the block (42), which is the 
output of the aforementioned set, is amplified by a factor NGAIN and multiplied by the signal 
produced by the envelope generator (43), whose time progression is illustrated in Fig. 14. The signal 
NOISE is the output of the generator (1 1) of Fig. 2. 

Fig. 12 describes the non-linear block "RATE LIMITER" (42) formed by the sums RLS1, RLS2, the 
limiter (44) and the unit delay element (45). The value memorised in the delay (45) is subtracted by 
means of the adder RLS1 from the input signal "IN"; the result is then limited to values within the 
interval ±RATE (being RATE the sequence generated by the network illustrated in Fig. 10), and 
finally summed again to the current delay value (45) at the node RLS2. The result "OUT" is 
memorised in the delay element (45) for a successive cycle. Fig. 13 shows the time progression of the 
envelope generated by the block (39): upon a "note on" event, starting from the level NBL0, the level 
NBL1 is reached in a time NBT, indefinitely sustained over time, also after the corresponding "note 
off" event. Fig. 14 shows the time progression of the sequence generated by the block (43): upon a 
"note on" event, the signal starts from the value NL0, reaches the value NL1 over a time NT1 and the 
level NL2 over a time NT2 sustained until the successive "note off" event. Upon this event the signal 
reaches the value zero in a time NT3. 

With reference to Fig. 1 1 , the non-linear block "RATE LIMITER" (42) can be replaced with a linear 
filter, whose gain has a progression described by the same sequence RATE generated by the 
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architecture of Fig. 10, so that the structure "NOISE BOX" of Fig. 1 1 is a linear time-variant filter. 
With reference to Fig. 2, the outputs of the generators (9) and (1 1) are the inputs of the resonator (12) 
illustrated with details in Fig. 15. The functional blocks of the network (12) form a cycle of operations, 
along which a sequence of samples propagates for a potentially infinite time. The two contributions of 
the two generators (9) and (1 1) are added to this sequence, instant by instant in the sum nodes (46) and 
(48) nodes, respectively, to sustain the energy of the computed sequence. The structure of Figure 15 is 
the translation into a mathematical model of the resonant part of the organ flue pipe, defined as 
"pipework". In particular, the low-pass filter (47) emulates the dissipation of acoustic energy, with 
variable intensities in function of the frequency; the high-pass filter (49) attenuates all the frequency 
components lower than the fundamental frequency; by means of the product (51), the envelope 
generator (50) produces a signal that represents the time progression of the loop gain of the resonant 
system; the filter (52) alters the sequence phase, leaving its module unchanged; the factor TFBK (53) 
depends on the type of acoustic termination at the top of pipework; finally, the delay line BDELAY 
(54) considers the time needed by an acoustic pressure wave to cover the pipework from the base to the 
top and vice versa. The time progression of the signal produced by the envelope generator (50) is 
traced in Fig. 16: likewise the envelope of Fig. 13, upon a "note on" event, the signal passes from a 
value FBL0 to a value FBL1 in a time FBT, and then remains constant. The output sequence (13) is the 
signal emitted by the mathematical model of Fig. 2 as a whole, that is to say the time representation of 
the sound emission of the organ flue pipes. 

The description continues with the original innovative characteristics of the audio digital synthesis 
technique of the sound of flue pipes. 

The literature on the generation of sounds of instruments with continuous sound emission, among 
which aerophone instruments, by means of the physical modeling technique, proposes solutions based 
on a mutual interaction between a non-linear active part, normally defined as excitation (55), and a 
linear passive part, defined as resonator (56), according to the scheme of Fig. 17. In the case of 
aerophone instruments, the energy contributed to the system is in the form or sound pressure and the 
signal produced is the progression of the sound pressure wave irradiated by one or more suitable points 
of the resonator. The waveform p(t) is the progression of the air pressure that the performer (or the 
bellows, in the case of a church organ) exercises on the instrument mouthpiece. According to this 
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progression and to the progression of the pressure w(t) in a suitable point inside the resonator, an 
oscillating acoustic pressure e(t) injected in the resonator is generated. Once the sustain phase has been 
reached, the pressure e(t) has the same fundamental frequency as the pressure w(t). Being linear 
(except for very special operation modes), the resonator can be described with an impulse response 
r(t), which generates the return signal w(t) and an impulse response h(t), which generates the output 
signal y(t). The latter is the time progression of the sound emission of the instrument. Being it a 
numerical simulation performed in the time domain instead of the frequency domain, the fundamental 
frequency of the oscillation on which the system stabilises, once the sustain phase has been reached, is 
extremely difficult to predict mathematically. This depends on the fact that the frequency depends on 
the time progression of the forcing signal e(t), and not only on the frequency values in which the 
amplitude spectrum of the impulse response of the resonator has the relative maximum values. In fact, 
any type of harmonic oscillator (electronic, mechanical, etc.) has this characteristic. With regard to 
wind instruments (including organ pipes), it is sufficient, for example, to increase the sound pressure 
to obtain an increase of the fundamental frequency of the acoustic wave, in addition to an intensity 
increase, although the characteristics of the resonant part remain unchanged. 

Another inevitable characteristic of the oscillating systems illustrated in Fig. 17 is the unpredictability 
of the phase of the generated signal, once the sustain phase has been reached. Since the waveform p(t) 
used to stimulate the system is partially chaotic, and in any case it does not contain any information 
about the phase of the stationary wave sustained by the resonator, the attack transient of the signal y(t) 
is always and unpredictably different. Therefore, although the waveform has always the same periodic 
time progression in sustain conditions, it is impossible to determine the evolutions that bring the 
system towards this progression. In quantitative terms, it is impossible to determine the phase of the 
fundamental frequency of any signal processed inside the stable oscillating system of Fig. 17, taking 
the instant when the stimulus p(t) starts as time origin. Together with the difficulty encountered in 
determining the fundamental frequency a priori, this is unacceptable in the field of high-polyphony 
electronic musical instruments, such as church organs. 

The synthesis system used in the present invention derives directly from the synthesis in the time 
domain described in general and is characterised by the total autonomy of the excitation signal from 
the signal produced by the resonator. In fact, the main harmonic sequence (10) extemporarily 
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generated by the block (9) of Fig. 2 is the imitation, as faithful as possible, of the signal e(t) of the 
system of Fig. 17 (assuming that the latter is a good mathematical model of the flue pipe of a church 
organ), with the substantial difference that the fundamental frequency and the phase of this sequence, 
and consequently of the sequence produced by the system (13) as a whole, are perfectly determined a 
priori. 

The preparation of the numerical parameters of any oscillating system, as generically illustrated in Fig. 
17, requires special sensitivity and skill, apart from the perfect knowledge of its mathematical model. 
This means that the good operation of the system may be impaired, and the system may become 
unstable or even inharmonious, if only one of the parameters has a value not included in a proper 
range. Moreover, different operational modes of the oscillator can be obtained only by acting 
simultaneously and with special attention on a plurality of parameters, with the risk of making the time 
evolution of one or more signals in transit along the functional blocks of the system uncontrollable. 
This makes the search for multiple sounds produced by this type of synthesis slow and difficult. On the 
contrary, a system with no feedback between resonator and excitation, such as the system shown in 
Fig. 2, enables to modify the numerical parameters of the three functional blocks (9), (1 1), (12) in a 
completely independent way, without impairing the good operation of the system as a whole. This 
allows obtaining a larger variety of sounds than the one obtained by means of a feedback loop system 
with equal complexity. 

The system of Fig. 3 shows a sequence of operations performed on the signal produced by the 
sinusoidal oscillator (14). The type and order of the operations are only one of the possible realisations 
used to generate a waveform sufficiently rich in harmonic components and provided with a suitable 
time evolution. In any case, some of the functional blocks of the system, such as the delay (24) and the 
non-linear function (26), derive from mathematical models of wind instruments known in the 
literature, without the need of using them. The originality of the system mainly consists in the 
adaptation of an ordinary oscillator with status variable to non stationary operational conditions, by 
developing the functional blocks (30) and (32) of Fig. 4, in order to make the oscillator robust to the 
variations of the parameter F 2 of the block (29). 

With reference to Figures 5 and 6, the originality derives from the development of the generator (34) to 
obtain pleasing random frequency variations in real time. Assuming the factor (33) as constant, that is 
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to say assuming the absence of the low frequency oscillation of the sequence 8f, the latter assumes a 
new random value at every period of the sinusoidal sequence VAR1 . The result is a statistic uniformly 
distributed variation of the wave period, in terms of probability density function. The variation is 
perceived as a pleasant irregularity in the sound emission. Otherwise, if 8f assumes a new random 
value at every sampling instant, the length of every wave period will be described by a variable formed 
by the sum of N aleatory contributions, each of them provided with uniformly distributed probability 
density (N is the number of samples per period). In view of the Central Limit Theorem, the higher is N, 
the more the probability density function of this variable approaches a Gaussian function. The 
frequency variation would be very irregular, since high frequency deviations would be obtained much 
more rarely than small deviations from the nominal frequency. This would be very unpleasant, since 
wave periods with very different length from the nominal length could be generated and perceived as 
sudden malfunctions of the generation model. 

The generator of the aleatory component (11) of Fig. 2 is completely original, and the embodiment of 
Fig. 10, 1 1 and 12 derives from the analysis of samples of sounds emitted by a large variety of organ 
flue pipes, and from some hypotheses on their operation physics. In particular, by analysing the 
spectrogram of the individual wave periods of a sample and using a much finer time resolution than a 
wave period, it can be noted that a large percentage of sound energy concentrates in a time interval 
much shorter than the period, always situated in the same position along the wave period. Such sound 
energy covers a frequency interval considerably higher than the interval covered as an average by a 
plurality of periods. Therefore, the characteristics of the spectrogram of the stationary part of the 
sound of flue pipes are similar to the spectrogram of a train of equidistant impulses, with the energy of 
the individual period concentrated in each impulse. These considerations justify the architecture 
illustrated in Fig. 10: the sequence RATE is obtained through a series of elementary deterministic 
operations performed on the sinusoid (16). Once the sustain phase has been reached, the sequence 
RATE assumes a qualitatively impulsive progression, of which Fig. 18 shows one example, where To 
is the period of the sinusoid (16). Regardless of the method used to obtain a pseudoimpulsive 
sequence, the sequence is conceptually one of the inventive foundations of the generator (11). 
The structure illustrated in Fig. 1 1 is formed by the four delay lines NBDL1, NBDL2, NBDL3, and 
NBDL4. Together with the sums NBS1, NBS2, and the product NCGAIN, the first three delay lines 
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form a FIR filter. The output of this filter (that is to say the sum NBS2) is processed by the non-linear 
element (42) and then, after being multiplied by NBFBK and after passing through the fourth delay 
line and the sum NBS3, injected again in the aforementioned filter. If it weren't for the element (42), 
the structure "NOISE BOX" would be a linear filter, whose spectrum would have a voluntarily 
inharmonious progression, with a plurality of resonance peaks distributed in a non-deterministic way, 
depending on the length of the delay lines and the two independent parameters NCGAIN and NBFBK. 
These four quantities are dimensioned in order to imitate the frequency response of a resonator with 
irregular geometry, such as the portion of space of the organ pipe immediately inside the mouth. 
Because of the periodic oscillation of the sequence RATE, the element (42) causes a continuous 
periodic variation over time of the "gain" (not strictly definable as such, since the "RATE LIMITER" 
is a non-linear block) of the entire "NOISE BOX". In particular, with reference to Fig. 1 8, when the 
sequence RATE assumes the minimum value, the non-linear distortions caused by the block (42) 
imply energy losses that heavily reduce the resonance effects of the "NOISE BOX". Vice versa, 
during the (much shorter) instants in which the sequence RATE assumes relatively high values, the 
resonant effect of the "NOISE BOX" emerges and the intensity of the aleatory component increases. 
It can be noted that during the attack transient, because of the envelope generator (39), whose 
progression is shown in Fig. 13, the sequence RATE assumes higher values than during the sustain 
phase; this increases the resonance of the "NOISE BOX" during the first instants of synthesis, in order 
to simulate the acoustic phenomena defined as chiff, cough, etc. produced by the flue pipes if the valve 
that regulates the passage of air from the bellows to the foot is opened rapidly. 
The non linear block (42) is formed by the two adders RLS1, RLS2, the limiter (44) and the unit delay 
element (45). At every sampling instant, the difference between the previous output value and the 
current input value is first limited in width to values within the interval +RATE and summed again to 
the previous output value, thus obtaining the current output value. The output sequence "follows" the 
input sequence, maintaining an inclination limited according to the value RATE. For mere illustrative 
purposes, Fig. 19 shows a chart of an input sequence (continuous line) and an output sequence (dotted 
line). In the instant to the inclination of the sequence IN exceeds the value RATE/sample, therefore the 
sequence OUT separates until it re-joins at point ti, after which the sequence IN remains constant. In 
the instant t 2 the excessive inclination of the sequence IN causes the immediate separation of the 



13 



Serial No. 10/549,388 
Docket No. 024783.003 

sequence OUT up to the re-conjunction point U. With respect to a linear filter, the advantage of the 
"RATE LIMITER" is the elimination of possible discontinuities of the aleatory sequence, while still 
maintaining a sufficient bandwidth, being such discontinuities extremely unpleasant for the human 
hearing. This aspect represents the originality of the "RATE LIMITER". 

The non linear block (42) can be replaced with any functional block whose effect on the structure 
"NOISE BOX" of Fig. 1 1 is the quantitative resonance variation generated by the structure, according 
to a periodic progression. 

As regards the linear resonator (12), the physical considerations that involve the choice of the 
functional blocks of Fig. 15 are described herein. The resonant part of an organ pipe, defined as 
pipework, can be mathematically described, in the most elementary way, with a "comb" filter 1/(1 - 
FBK-z ~ N ), in which the feedback coefficient FBK is related to the loop gain of the filter and the 
parameter N is inversely proportional to the first resonance frequency of the same. The more complex 
resonator of Fig. 15 derives from this base, which is very used in the field of audio digital processing. 
Among the elements of the resonator, the function of the delay line (54) appears evident. The response 
in module of the low-pass filter (47) is designed so as to consider the different energy losses suffered 
by the various harmonic components during their transit along the pipework, while the high-pass filter 
(49), whose cut-off frequency is lower than the fundamental frequency of the resonator, completely 
eliminates the continuous component of the stationary wave, to take into account the fact that the 
average acoustic pressure inside a pipework is approximately equal to the external pressure. Because 
of the envelope generator (50), during the first operation phase of the resonator, the loop gain of the 
system is moderately overabundant with respect to the value once the sustain phase has been reached, 
in order to obtain a faster initial energy accumulation in the resonator, that is to say a faster attack 
transient of the generated sound. The sign of the factor TFBK (53) is especially important: a positive 
sign for a pipework open at the mouth and on top, and a negative sign for a pipework open at the mouth 
and closed on top. This derives from the physics of the reflection of an acoustic pressure wave in 
correspondence of the pipework terminations. This physical law also justifies the use of the all-pass 
filter (52), the most important element of the resonator from the conceptual point of view. If, on one 
side, the mono-dimensional model of the pipework is sufficiently accurate to justify the use of an 
individual delay line to simulate the longitudinal propagation of an acoustic wave in the pipework, the 
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approximation becomes unacceptable in the wave reflection in correspondence of a structural 
discontinuity characterized by non-negligible transversal dimensions, such as the top of the pipework. 
The all-pass filter (52) modifies the total phase delay of the closed cycle formed by the elements (46) 
... (54) in a selective way with respect to the frequency, in order to make the resonance of the linear 
resonator (12) realistically inharmonious. The same filter is optionally used to modify the value of the 
first resonance frequency of the pipework in real time through controlled variations of its coefficients 
upon a "note off event, to simulate the phenomenon of the moderate reduction of the fundamental 
frequency of the sound of small flue pipes when the air inlet valve closes. 
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